The polydisperse cell model: 
Non-linear screening and charge renormalization in colloidal mixtures 
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We propose a model for the calculation of renormalized charges and osmotic properties of mixtures 
of highly charged colloidal particles. The model is a generalization of the cell model and the notion 
of charge renormalization as introduced by Alexander and his collaborators (J. Chem. Phys. 80, 
5776 (1984)). The total solution is partitioned into as many different cells as components in the 
mixture. The radii of these cells are determined self-consistently for a given set of parameters from 
the solution of the non-linear Poisson-Boltzmann equation with appropriate boundary conditions. 
This generalizes Alexanders's model where the (unique) Wigner-Seitz cell radius is fixed solely by 
the colloids packing fraction. We illustrate the technique by considering a binary mixture of colloids 
with the same sign of charge. The present model can be used to calculate thermodynamic properties 
of highly charged colloidal mixtures at the level of linear theories, while taking the effect of non-linear 
screening into account. 
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I. INTRODUCTION 



Typical charge-stabilized colloidal suspensions consist 
of mesoscopic particles and ions, immersed in a sol- 
vent which to a good approximation can be consid- 
ered as a continuum medium, characterized by a di- 
electric permittivity e. The microscopic ions form a 
diffuse layer around the colloidal surfaces, with a typ- 
ical thickness of the order of the Debye length of the 
electrolyte K , and with a net charge that compen- 
sates the charge of the colloidal surface. As a conse- 
quence, the effective colloidal interactions are screened: 
instead of the bare Coulomb interaction Z 2 e 2 /er between 
two colloidal spheres of radius a and charges Ze a dis- 
tance r apart in the solvent, standard linear Poisson- 
Boltzmann theory predicts a screened coulomb interac- 
tion Z 2 e 2 exp(2/ca)(l + na)~ 2 exp[— Kr)/er between the 
colloids [1, 2], However, it is well-known that this rela- 
tively simple picture needs modification in the case of 
highly charged colloids, where the strong electrostatic 
coupling between the colloids and their counterions in- 
duces the accumulation of the latter close to the colloid- 
surface. This condensation phenomenon cannot be de- 
scribed by a linearized Poisson-Boltzmann approach, and 
requires the full nonlinear Poisson-Boltzmann theory for 
a proper description [3], However, the simplicity of the 
linear theories can be retained by the use of the concept 
of charge renormalization, which considers each colloid 
and its condensed counterion shell as a single entity, car- 
rying a net effective charge which is usually smaller than 
the bare (structural) charge Z of the colloid. This so- 
called renormalized charge, Z len , characterizes the long- 
distance behavior of the screened electric potential cre- 
ated by the colloid and its condensed microion polar- 
ization cloud, and the effective colloidal pair interaction 



takes the linear-screening form but with Z replaced by 

Z ren 

The concept of charge renormalization, and also its 
relatives such as dressed- ion theory [4], is very well- 
developed by now for monodisperse colloidal suspensions 
[3, 5-7]. The traditional theory by Alexander et al. is 
based on a cell model, in which a suspension of N col- 
loidal spheres in an electrolyte of volume V is approxi- 
mated by N identical spherical cells of volume V/N, each 
containing a single colloid in the center and co- and coun- 
terions such that each cell is neutral. The simple spher- 
ically symmetric geometry allows for a straightforward 
numerical solution of the nonlinear Poisson-Boltzmann 
equation in the cell, from which Z lcn follows directly 
[3, 5-7]. Although such an approach based on cells 
is mostly justified for colloidal crystals, the concept of 
charge renormalization has proved to be very useful for 
fluid-state suspensions, e.g. in describing the structure 
factor [8] or the electrophoretic mobility [9], 

Much less is known about charge renormalization in 
the case of colloidal mixtures or in polydisperse colloidal 
suspensions. In recent years such mixtures of charged 
colloids have generated a lot of interest, for instance be- 
cause of their ability to form a plethora of crystal phases 
[10, 11] or fluid-fluid demixing phenomena [12]. More- 
over, in the presence of external fields colloidal mixtures 
also show a richer behavior than that of monodisperse 
systems, e.g. oppositely charged colloids can exhibit lane 
formation in an electric field [13-15], while sedimentation 
in the Earth's gravity field can generate layers of colloids 
with equal mass-per-charge [16] or a "colloidal Brazil nut 
effect" in which a layer of heavier colloids floats on top of 
a layer of lighter ones [17]. Given the success of the cell 
model and of the concept of charge renormalization for 
the calculation of thermodynamic properties of monodis- 
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perse highly charged colloidal suspensions, as well as the 
recent interest and advances in the synthesization of col- 
loidal mixtures, it is desirable to extend the cell model 
to the polydisperse domain. The aim of this paper is 
to introduce such a generalized cell model and to study 
some of its features, in particular the dependence on the 
composition of the mixture — this plays of course no role 
in monodisperse suspensions. 

The outline of this paper is as follows. In section 2 we 
introduce the model. In section 3 we discuss the partic- 
ular case of a binary mixture. We calculate saturation 
curves for the renormalized charges, osmotic properties, 
and the radii of the cells as a function of composition, 
screening length and total colloid density. Discussion and 
some general remarks are presented in section 4, followed 
by conclusions in section 5. 



II. MODEL 

Let us begin by considering two species of spherical 
colloids with charges Q a = Z a e, Q b = Z b e and radius 
a and b, respectively, where e is the proton charge and 
where we assume that the colloids possess charges of the 
same sign, i.e., Z a Z b > 0. Let N a and N b denote the total 
number of colloids of type a and those of type b respec- 
tively. The colloidal mixture occupies a volume V and 
the colloids are immersed in a medium that is in osmotic 
equilibrium with a salt reservoir with screening length 
= (8ttp s \b)~ 1 / 2 where 2p s is the salt concentration 
in the reservoir. As usual we have defined the Bjerrum 
length Xb = f3e 2 /e, with s the dielectric constant of the 
solvent and /3 = l/(fcsT), with T the temperature and 
ks the Boltzmann constant. The average number densi- 
ties of colloids are defined by p a = N a /V and p b = N b /V 
and their respective volume fractions are r/ a = Ana p a /3 
and rj b = 4irb 3 p b /3. The composition of the mixture is 
characterized by the molar fraction x = N a /(N a + N b ) 
which in the case of equisized colloids can be written as 
x = r) a /rj with r] = r\ a + rj b the total packing fraction of 
the suspension. 

Following the usual cell model approximation we con- 
sider the Wigner-Seitz cell of each colloidal particle as 
a sphere concentric with the spherical colloid. However, 
due to the difference in size and charge of each species, 
there is no a priori reason to assume that the cells cor- 
responding to each species are of equal size. Therefore, 
we suppose that the Wigner-Seitz cells of the colloids of 
type a have a radius R a and those of the colloids of type 
b have radius R b , which are in principle different and to 
this point unknown. 

Let <p a {r) = (3eip a (r) be the reduced electric potential 
inside a cell of type a, with ip a (j) the electric potential 
at a radial distance r from the center of the colloid. Ac- 
cordingly, we define 4* b (r) = (3eip b (r) the reduced electric 
potential inside a cell of type b. In the mean-field ap- 
proximation, the electric potential inside each cell is the 



solution of the Poisson-Boltzmann equation 

A<f) a (r) = k 2 sinh</> a (r) 
A(j) b (r) = k 2 sinh (j) b (r) 

together with the boundary conditions 



(la) 
(lb) 



dr 



d<j)g 

dr 



= -Z a \ B /a; r 
= 0; 



r=R a 



dr 
d(f> b 
dr 



= -Z b X B /b (2) 



r—b 



r=R b 



(3) 



which are obtained from direct application of Gauss' law 
at the colloid surface (Eq. (2)), and at the boundary 
of the cell (Eq. (3)). The latter equation accounts for 
the electroneutrality of each cell. These equations com- 
pletely determine the electrostatic potential inside each 
cell. However, the radius of each cell R a and R b remain 
undetermined, so far. These are determined by the fol- 
lowing relations: 



R 3 R 3 
a 6 b 6 



(4) 



which expresses the fact that the Wigner-Seitz cells fill 
the whole accessible volume V, and 



4>ai.Ra) = <f>b(Rb) , 



(5) 



which is a consequence of the continuity of the electric 
potential, and the fact that cells of different species can 
be in contact with each other. The system of differen- 
tial and algebraic equations (1), (2), (3), (4) and (5) is a 
complete set to determine the unknown parameters and 
functions R a , R b , (pair) and </>b(r). The numerical so- 
lution of this system of equations can be obtained with 
a generalization of the elegant algorithm for a monodis- 
perse colloidal suspension proposed in Ref. [5]. 

The generalization to a mixture of M > 2 colloidal 
species with radii r i; packing fractions rji and charges 
Zi , is performed by introducing M cells of radii Ri (i = 
1, 2, M), in which the Poisson-Boltzmann equation (1) 
with boundary conditions (2) and (3) is to be solved. To 
fix the values of the M radii of the cells, the equivalent 
of Eq. (4) is 



M 7?3 



(6) 



Additionally, we have M — 1 independent equations of 
continuity for the electric potential at the cell boundaries 
of the form 



i(Ri) = MRi)- 



(7) 



for i = 2, M. This completes the system of equations 
for the electric potentials and the radii of the cells. 
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A. Renormalized charge 

The renormalized charge, defined using Alexander pre- 
scription [3], is obtained by comparing the nonlinear so- 
lution of the Poisson-Boltzmann equation with its lin- 
earized solution, where the linearization is performed 
with respect to the value of the potential at the bound- 
ary of the cell. Following Ref. [5], and restricting the 
analysis to a binary mixture for the sake of simplicity, let 
us define 4>i(r) — (f>i(r) — (f>o, for each cell i = a,b, with 
4>a = 4>a{Ra) = (pbiRb)- Supposing that <j>i is small, it 
satisfies the linearized Poisson-Boltzmann equation 

A4>i(r) =4 B (^(r)+ 70 ) (8) 



with 



k 2 cosh tfio and 70 = tanh < 



ary conditions 4>i{Ri) = and (^[{Ri) 
of the linearized problem is 



, and the bound- 
0. The solution 



<f>i(r) = 7o 



- 1 



3 K PB (-Ri 



+- 



e^' r - fl -)( K p B ^ + l) 
2np B r 



2k pb t 



(9) 



The linear potential <fii(r) + 4>o, the behavior of which 
approximates the nonlinear solution 4>i{r) near the cell 
boundary, is created by an effective or renormalized 
charge Z\ cn = r^ip'^n)/ Xb where is the colloid radius 
(r a = a, rb = b for i = a, b respectively). Explicitly, 

L ' ^(«;| B i?irj - l)sinh(K PB (-Ri - n)) 

n) cosh(Kp B (Ri - n)) . (10) 



Kp B r 

+ K PB (Ri 



As mentioned before, a convenient numerical algorithm 
to compute the renormalized charges can be formulated 
on the basis of the method introduced in Ref. [5]: we 
impose a value </>o for common value of the electric po- 
tential at the cell boundaries. Solving numerically the 
system of equations (l)-(5) allows us to determine the 
radius of the cells R a and Rb and the bare charges Z a 
and Zt corresponding to the value </> of the potential 
at the boundaries of the cells. Finally, using Eq. (10) 
we compute the corresponding renormalized charges Z a en 
and Z v b en . The osmotic properties of the suspension can 
also be determined by this procedure as discussed below. 



B. Thermodynamic properties 

In the mean field approximation, the grand potential 
per cell (w a or LOb) is given by (see e.g. Ref [18]) 

fai = 7^T~ / |V<Mr)| 2 dr 



8tt\ b 
+ 



I 



1) 



a=± 

f3no[pt(r)+pr(r)-2p s ]dr, 



dr 



(11) 



where the ionic density profiles inside each cell are given 

by 



pf(r)=p s e^. 



(12) 



The integration is performed in the cells excluding the 
volume occupied by the colloid, i.e., rj < r < Ri- 

Using Gauss' law, Eqs (2) and (3), this can be written 

as 



Zi 



(13) 



rRi 

+ 47T p s I dr r 2 [(f>i(r) smh 4>i{r) — 2 cosh (f>i(r)] 

J Ti 



The total grand potential of the system is 

n = N a cu a + N b Lo b + K d + n d 



(14) 



where f3Ff = A^hi^A?) - 1] accounts for the colloid 
entropy, with Aj and iVj the thermal length of colloids of 
species i and their number respectively The electrostatic 
contribution to the pressure p can be obtained deriving 
ft with respect to the volume V. Introducing the os- 
motic pressure Yl= p — p res with p res the pressure in the 
electrolyte reservoir, the final result can be expressed in 
terms of the potential at the boundary of the cells 0o in 
the form 



^E = l + ^£ ( C o S h O -1) 
P P 



(15) 



with p = p a + Pb the total number density of colloids. 
The l.h.s. of Eq (15) defines the osmotic compressibility 
factor which equals one in the dilute limit. In the next 
section we apply the model formulated above to study 
charge renormalization and osmotic properties in a col- 
loidal binary mixture. 



III. COLLOIDAL BINARY MIXTURE 

Let us consider a colloidal mixture consisting of two 
species of colloids. In order to reduce the number of pa- 
rameters, we will assume that the colloids have the same 
radius, i.e., a = b and we will express all lengths in units 
of this quantity. In addition, the Bjerrum length will be 
fixed to \b I a = 0.0022 which is a typical value for a 
suspension of colloidal particles in water. The remaining 
free parameters are (i) the bare charge of species a: Z a , 
(ii) the bare charge of species 6: Zb (iii) the total pack- 
ing fraction: 77, (iv) the composition of the mixture: x, 
and (v) the concentration of electrolyte as given by the 
parameter na, where k _1 is the screening length of the 
electrolyte reservoir. Notice that the screening parame- 
ter k is different from kpb as defined after equation (8), 
the later corresponding to the electrolyte concentration 
at the cell boundary. 



4 



A. Renormalized charges 



Let us start our analysis by considering the behavior 
of the renormalized charges of the binary mixture as a 
function of the structural charges. Fig. 1 shows the renor- 
malized charge for each species when the bare charge of 
one of the species is fixed while that of the other species 
is changed. In Fig. la the structural charge of species b 
is fixed to Z^Xs/a = 5 while Z a varies. In Fig. lb we 
fixed Z a Xs/b = 10 and vary the bare charge of species 
b. Two different compositions are considered, x = 0.05 
and x — 0.95, and two screening lengths na — 0.1 and 
na = 1. The remaining parameters are fixed to r\ = 0.1, 
Xb = 0.72nm and a = b = 326nm. In Fig. 2 a similar 
study is performed, this time for a larger total packing 
fraction of the mixture rj = 0.3. A first observation fol- 
lowing Figs. 1 and 2 is that the dependence of Z\ cn is 
similar to the behavior found in monodispcrse systems: 
as the bare charges Zi increases, the corresponding renor- 
malized charges Z\ en grows. The growth is only linear for 
small values of the varying structural charges. For larger 
values of ZiXs/a the increment becomes highly non- 
linear (compared to the dashed lines corresponding to 
Z[ e " = Zi) and saturates towards a value much smaller 
than the structural charge for ZiXs/a ^> 1. Notice that 
the saturation value of the charges depends not only on 
the total packing fraction but also on the composition 
of the mixture. The influence of the electrolyte concen- 
tration on the saturation curves, on the other hand, is 
relatively small, as compared to the effect produced by 
a drastic change in composition or total packing fraction 
of the mixture. In particular, in Fig. 2 the curves corre- 
sponding to na = 1 and na = 0.1 arc indistinguishable. 
For the packing fractions considered in these examples, 
reducing the electrolyte concentration manifest itself in 
a slight reduction of the saturation values in all cases. 

The concentration of colloidal particles in suspensions 
can be adjusted and measured in a wide range of pack- 
ing fractions, therefore, it is interesting to analyze to 
what extent the general behavior of the renormalized 
charges described above depends on this quantity. In 
Fig. 3 we consider the scaled renormalized charges cor- 
responding to a total packing fraction varying over sev- 
eral decades. The bare charges of species a and b are 
fixed by Z a As/a = 10 and Z^Xs/a = 5 respectively. 
The figure shows curves obtained for different electrolyte 
concentrations as fixed by the parameter na. Several 
values of a; = r] a /r) are represented by different types of 
lines. For na > 1 and r\ < 0.01 varying either the elec- 
trolyte concentration, the total packing fraction, or the 
composition of the mixture has a small effect on the val- 
ues of the renormalized charges. For r\ > 0.01 varying 
the composition of the mixture causes the charge renor- 
malization curves to split into several branches. This 
effect is also present for smaller total packing fractions 
in the extremely low electrolyte concentration curve cor- 
responding to na = 0.1. In the latter case, the values 
of the renormalized charges are much more sensitive to 
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FIG. 1: Renormalized charges Z„ en As/a and Z^ en As/a for 
a binary mixture with: (a) fixed bare charge i^As/a = 5 
and varying bare charge Z a Afi/a and (b) fixed bare charge 
ZaXs/a — 10 and varying bare charge Z 6 As/a for Xb = 
0.72nm and a = b = 326nm. The total packing fraction is 
fixed to the value r\ = 0.1. Different electrolyte concentra- 
tions and compositions are considered as pointed by the leg- 
ends. For the sake of comparizon, the dashed line shows the 
linear dependence corresponding to the case Z\ en — Zi. 

variations of the total packing fraction, which increases 
after reaching a minimum at about 77 = 0.015. At total 
packing fractions 77 > 0.2 all the curves show the same 
qualitative behavior and tend to merge into lines of sim- 
ilar slope. 

In Fig. 4 we show the effect of varying the composi- 
tion of the mixture on the scaled renormalized charge 
of species a. We observe that ZJ en As/a varies linearly 
with changing x between the corresponding values of a 
pure species— a system, i.e., x = 1 and a pure species— b 
system, i.e., x = 0. Similar curves are obtained for 
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FIG. 2: Renormalized charges for a binary mixture as given 
in Fig. 1. Here we consider a higher value of the total packing 
fraction, rj = 0.3. 

Z[ en \B/a. The slope of the renormalized charge curves 
is slightly larger for small electrolyte concentration and 
is practically zero for na — 2. For the larger packing 
fractions (triangles), the general tendency is that the 
renormalized charges increase with x. For smaller pack- 
ing fractions and low electrolyte concentration (squares 
in Fig. 4a), this tendency is reversed in correspondence 
with the behavior observed in Fig. 3a. 



B. Osmotic properties 

In Fig. 5 we show the compressibility factor (Eq. (15)) 
as a function of the total packing fraction, correspond- 
ing to the four different electrolyte concentrations con- 
sidered in Fig 3. Different types of lines illustrate the 



FIG. 3: Renormalized charges Z^Xb/o, and Z£ en As/a for 
a binary mixture with bare charges ZuXb/o, = 10 and 
Z b \B/a = 5 as a function of the total packing fraction 77. 
Different electrolyte concentrations as fixed by the parameter 
na are shown. The values of the remaining parameters are 
typical for colloidal suspensions, namely Xb = 0.72nm and 
a = b = 326nm. Different types of lines correspond to dis- 
tinct compositions of the mixture as indicated by the legends. 

effect of varying the compositions of the mixture. In gen- 
eral, adding particles of species with the larger charge in- 
creases the compressibility factor. This effect is particu- 
larly notorious at high packing fractions. At low packing 
fractions the electrostatic pressure reduces to the value 
of the reservoir pressure and the compressibility factor 
decays towards the ideal gas value, i.e., /JIT/ p =1. It 
is interesting to observe the strong effect of varying the 
electrolyte concentration: following the sequence from 
Fig. 5a to Fig. 5d one notices that the compressibility 
factor decreases by about two orders of magnitude for 
rj < 0.1. This diminishing in the compressibility factor is 
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FIG. 4: Renormalized charge Z„ e "As/a as a function of the 
composition x = r/ a /rj for different salt concentrations for the 
parameters b — a, ZaXs/a = 10, and Z^Xsja — 5. We con- 
sider different values of the total packing fraction, namely 
r\ = 10~ 4 , (circles), r\ = 10~ 3 , (squares), r\ — 1CF 2 (dia- 
monds), and rj = 1CF 1 (triangles). 

a consequence of the decreased range of the repulsions. 
Moreover, for low electrolyte concentrations the com- 
pressibility factor rapidly grows as a function of the total 
packing fraction. This effect is intensified by increasing 
the concentration of the highest charged species, species 
a in this particular case study. For higher electrolyte con- 
centrations, Fig. 5d, the rate of change of the compress- 
ibility factor is reduced and it assumes smaller values 
as comparison with Fig. 5a shows, i.e., the system be- 
comes ideal-like in the dilute regime, as a consequence of 
the screening of the strong repulsive electrostatic forces. 
This can be seen more clearly from the osmotic equa- 
tion of state shown in Fig. 6 for several electrolyte con- 
centrations. In Fig. 6a, corresponding to extremely low 
screening, the osmotic pressure assumes large values and 
the suspension is expected to be in crystal phase. After 
increasing the screening constant, the osmotic pressure 
assumes values close to zero which means that the excess 
pressure due to electrostatic interactions reduces to the 
reference pressure in the reservoir. 



FIG. 5: Compressibility factor as obtained from Eq. (15) 
as a function of the total packing fraction for different elec- 
trolyte concentrations. The line-type convention for different 
compositions the same as in Fig. 3. 



0.05 0.10 0.05 0. 
— I ' 7~^i I ' 1 ' 1 




FIG. 6: Osmotic pressure as a function of the total packing 
fraction for different electrolyte concentrations. The line-type 
convention is the same as in Fig. 3. 



C. Cell radii 

In Fig. 7 we plot the cell radii (expressed in units of 
the colloidal radius a) as a function of the charge of one 
of the species while the charge of the other species is 
fixed (see details in the caption). In general, as one of 
the radii increases the other radius decreases to satisfy 
the constraint of fixed total packing fraction as expected 
from Eq. (4). The form in which this occurs depends 



on the colloid charges and compositions as well as on the 
other physicochemical parameters in the model. A small 
total packing fraction allows for larger cell radii as we 
notice by comparing Figs. 7a and 7c on the one hand 
and Figs. 7b and 7d on the other hand. Notice also that 
in the latter case for small values of Z^Xs/a we have 
R a /Rb > 1- After the charge of species b is increased, 
the inequality quickly inverts i.e., the volume of the cells 
surrounding the now highly charged particles of species b 
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FIG. 7: Radii of the cells f? and Ri, surrounding particles 
of type a and b respectively expressed in units of the colloids 
radius a = b = 326nm. In (a) and (c) the bare charge of 
type-a particles varies while Z^Xb/cl = 5 is kept fixed. In 
(b) and (c) the bare charge of type-b particles varies while 
ZaXs/a = 10 remains unchanged. Two different total packing 
fractions are considered as pointed by the legends and the 
values As = 0.72nm and kci = 1 are used in all figures. The 
line-type convention is the same as that in Fig. 3. 



surpasses that of species a. Finally, it is worth to notice 
the fact that beyond some value of ZiXs/a the radii of 
the cells become almost charge independent, reflecting 
the fact that the renormalized charges have reached their 
saturation value. In the light of the previous case study, 
we now proceed to discuss in more detail some general 
features of the polydisperse cell model. 



IV. DISCUSSION 

Physically, the charge renormalization process results 
from the accumulation of ions in the vicinity of a highly 
charged object immersed in an electrolyte. Therefore, the 
relevant parameter to compute their interactions is not 
the bare charge but the effective (renormalized) charge 
that takes into account this non-linear screening effect. 
As shown in Fig. 1, the behavior of the renormalized 
charges for polydisperse systems at the level of PB the- 
ory is reminiscent of that of their monodisperse coun- 
terpart: the renormalized charges saturate to a finite 
value when the structural charge becomes sufficiently 
large. The saturation values can be orders of magnitude 
smaller than the structural charges. As a consequence of 
this non-linear screening effect, a very asymmetric col- 
loidal mixture behaves as a much more symmetric sys- 
tem as regards their electrostatic effects. For instance, 
from Fig. 2b one can notice that a colloidal binary mix- 



ture with charges ZaXs/a = 10 and Z^Xb/cl = 300 be- 
haves effectively as a mixture with ^ en As/a ~ 7.5 and 
13, so that the ratio Z a /Zb = 0.03 changes 
= 0.58 which corresponds to a much more 
charge symmetric effective mixture. On the other hand, 
for small values of the structural charges, Z\ en ~ Zi. 
As a consequence, expressions for the effective interac- 
tion potentials derived from linear theories, e.g., DLVO 
theory, are expected to be sufficient to describe the ther- 
modynamic properties of the suspension, at least in the 
regime where pairwise additivity is expected [19, 20]. 

The saturation values of the renormalized charges de- 
pend on the physicochemical parameters: the electrolyte 
concentration, total colloids density, the dielectric con- 
stant of the solvent and the composition of the mixture. 
Increasing the electrolyte concentration has the effect of 
enhancing the value of the renormalized charges. This 
effect induces two competing tendencies in the screened 
interaction potential. On one hand the spatial extent 
of the double-layer reduces, on the other hand the am- 
plitude of the interactions increases. This tendency is 
inverted in the case of weak screening and total colloids 
packing fraction smaller that about 5% (see Fig. 3). The 
balance of these two non-linear screening features is af- 
fected by the composition of the mixture and is expected 
to have consequences on, e.g., the phase diagrams of mix- 
tures of highly charged macroions computed by using 
Yukawa-like potentials, the omission of such non-linear 
screening effects, as accounted by Z\ er \ may lead to un- 
physical results as demonstrated by several authors (See 
e.g. [21, 22]). 

In the context of the cell model for monodisperse sys- 
tems, it is assumed that due to entropic reasons, the 
colloidal particles tend to occupy the maximum space 
possible inside the solvent. Therefore, it is feasible to 
factorize the suspension into spherical cells, each of them 
housing a colloidal particle in its center. The radius of 
such cells is then chosen in such a fashion that the cells 



fill all the suspension, namely R/a 



V 



_1 / 3 . In other 



words, the packing fraction alone determines the size of 
the cells. In contrast, in the polydispense cell model, the 
radii of the cells depend not only on the total packing 
fraction and the composition of the mixture according to 
Eq. (4), but also on the relative charges of the species 
through the boundary conditions expressed by Eqs. (3) 
and (5). In other words, the volumes of the cells are 
not only determined by fixing the total packing fraction, 
but adapt themselves in order to satisfy the boundary 
conditions which impose the cells to be electrically neu- 
tral, the potential to be continuous at the cells bound- 
aries and the cells to fill the whole volume accessible to 
the colloidal particles. The resulting cells are indepen- 
dent units regarding electrostatic interactions in virtue 
of their neutrality. Therefore, the grand potential of the 
entire suspension can be easily calculated as the sum of 
the grand potential of individual cells. This facilitates 
enormously the calculation of thermodynamic response 
functions, in particular, following the contact theorem, 
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(see e.g., Ref. [2]), the osmotic pressure can be obtained 
directly from the value of the electrostatic potential at 
the cells boundary as explained in Sec. 2. The value of 
the potential at the cell boundary depends on the size of 
the cell, which have a priori different radii depending on 
the parameters, and in particular, on the composition of 
the mixture. 

Experimentally, in spite of the recent advances in elec- 
trophoretic techniques [23] , it is difficult to precisely mea- 
sure the values of renormalized charges in colloidal sus- 
pensions. Therefore, it would presumably be intricate 
to collect experimental evidence of, e.g., the splitting 
of the charge renormalization curves in Fig. 3, even at 
the highest packing fractions. Osmotic properties are 
nevertheless relatively easy to probe experimentally e.g., 
in sedimentation-diffusion equilibrium [24], [25]. In that 
sense the present formalism provides a reasonably sim- 
ple and numerically robust scheme for the calculation of 
the osmotic equation of state and compressibility factors 
in a wide range of parameters. As illustrated in the case 
study of Sec. 2, these osmotic properties are no n- trivially 
dependent on the mixture composition, which illustrates 
the fact that the relatively small change in the values of 
the corresponding renormalized charges shown in Fig. 3, 
do not necessarily imply a mild change in the response 
functions: the former are obtained from changes in the 
slope of the potential at the colloids surface as given by 
the solution of the linearized PB equation (see Eq. 10), 
whereas the latter depend on the value of the potential 
at the cells boundary (Eq. 15). The screening param- 
eter, however, plays a much more important role in the 
osmotic properties than the relative composition of the 
mixture, and may cause, for instance, the osmotic pres- 
sure to vary over several orders of magnitude. 

It is worth to mention the fact that in the polydisperse 
cell model we restricted the attention to colloids of the 
same sign of charge. The reason is that one can impose 
on a physical ground that each cell surrounding a colloid 
is charge neutral i.e., = for i = 1, M. It is not 

obvious that such a condition can be imposed for oppo- 
sitely charged colloids, as the potential may change sign 
in between the colloids while the electric field does not 
vanish anywhere. The attractive or repulsive nature of 
the interactions between the cells is not obvious to esti- 
mate by analyzing the sign of the potential, since e.g., 
potentials opposite in sign can still lead to colloidal re- 
pulsions [26]. A second limitation related to the cells 
electoneutrality, which is also present in Alexander's for- 
mulation, is that as a consequence of the restriction of the 
theoretical analysis to a single cell per colloidal compo- 
nent, all information on correlations between the colloids 
is lost. This may lead to unphysical ideal-gas- like values 
for the osmotic pressure in highly screened suspensions 
[27]. Care has to be exercised in drawing conclusion from 
cell models in this regime. 

An important feature of charge renormalization is the 
fact that increasing the electrolyte concentration leads to 
higher values of the renormalized charges simultaneously 



reducing the extension of the double-layer [6, 7]. These 
two non-linear screening features are neglected in linear 
theories, such as the DLVO model, where the (structural) 
charge is in principle independent of the screening pa- 
rameter. These artifacts of the linearization procedure 
can be cured by using the renormalized values of the 
charges within the linear schemes, which can in principle 
be achieved for example by using the analytical solution 
of the Ornstcin-Zernikc equation for a mixture interact- 
ing via a DLVO-like potential within Mean Spherical Ap- 
proximation, by replacing Zi — > Z\ m . A second alterna- 
tive is provided by the variational treatment of DLVO- 
like potentials based on the Gibbs-Bogoliubov inequality, 
which uses the Laplace transforms of the pair correla- 
tion functions to estimate the free energy of the system 
though a numerical optimization problem (see e.g., Refs. 
[28, 29]). This alternative is extremely suitable for the 
calculation of phase diagrams by e.g., using Maxwell con- 
structions, providing a simple alternative for the study of 
phase behavior of highly charged colloidal mixtures while 
fully taking into account non-linear effects. 



V. CONCLUSION 

In this paper we proposed a simple model which al- 
lows to study the thermodynamic properties of mixtures 
of highly charged macroions suspended in an electrolyte. 
The model is a generalization of the well-know cell model 
for colloids, and the notion of charge renormalization as 
introduced by Alexander and collaborators [3]. In con- 
trast to the one-component cell model, where the radius 
of the cells is fixed by the packing fraction of colloids, in 
the polydisperse cell model the radii of the cells depend 
on the physicochemical parameters through the bound- 
ary conditions utilized to solve the non-linear Poisson- 
Boltzmann equation in the cell geometry. The radii of 
the cells play the role of Lagrange multipliers that en- 
force the constraint that the cells fill the entire volume 
occupied by the colloidal mixture and have to be deter- 
mined within the numerical solution of the non-linear dif- 
ferential equations. We presented a method to perform 
such a calculation based on the algorithm of Ref. [5], 
which also allows to determine thermodynamic response 
functions and renormalized charges within a clear and 
efficient numerical scheme. We presented a detailed ex- 
ample of such an analysis in the case of a colloidal binary 
mixture, where we studied the effect that several physico- 
chemical parameters have on the values of the saturation 
curves for the renormalized charges and osmotic proper- 
ties of the suspension. In general, the ion condensation 
phenomenon that leads to charge renormalization has the 
effect of making the mixture more "electrostatically sym- 
metric" , in the sense that the ratio of the renormalized 
charges is much closer to one than that of the structural 
charges. Increasing the concentration of the species with 
the greatest charge keeping the total packing fraction 
fixed, raises the osmotic pressure. An important limi- 
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tation of the present model, which is actually shared by 
all Poisson-Boltzmann-based models, is the fact that ion- 
ion correlations are neglected. This oversimplification of 
the problem may lead not only to quantitative but also 
qualitative modifications to the saturation curves: the 
renormalized charge may exhibit a maximum followed 
by a decrease of Z\ en with Zi increasing (see e.g., Ref. 
[30-32]). In some cases this can lead to a weakening of 
the repulsive effective interactions as calculated on the 
basis of the renormalized charges, which may eventually 
turn electrostatic repulsions into attractions. In the case 
of 1:1 electrolytes, such ion- ion correlations are relatively 



unimportant and the present Poisson-Boltzmann model 
is expected to provide a reasonably good description of 
highly charged colloidal mixtures. 
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